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O ■ ABSTRACT 



Neither standard model SEDs nor truncated standard model SEDs fit ob- 
served spectra of QU Carinae with acceptable accuracy over the range 900A to 
3000A. Non-standard model SEDs fit the observation set accurately. The non- 
standard accretion disk models have a hot region extending from the white dwarf 
to R = 1.36i?wd, a narrow intermediate temperature annulus, and an isothermal 
remainder to the tidal cutoff boundary. The models include a range of M values 
between 1.0 x 10" ''Mq yr"^ and l.OxlO~®M0 yr"^ and limiting values of M„d 



^ Based on observations made with the NASA/ESA Hubble Space Telescope, obtained at the Space 
Telescope Science Institute, which is operated by the Association of Universities for Research in Astronomy, 
Inc. under NASA contract NAS5-26555, and the NASA-CNES-CSA Far Ultraviolet Explorer, which is 
operated for NASA by the Johns Hopkins University under NASA contract NAS5-32985 
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between O.6M0 and 1.2Mq. A solution with M^d = I.2M0 is consistent with 
an empirical mass-period relation. The set of models agree on a limited range of 
possible isothermal region T^a values between 14,000K and 18,000K. The model- 
to-model residuals are so similar that it is not possible to choose a best model. 
The Hipparcos distance, 610 pc, is representative of the model results. The 
orbital inclination is between 40° and 60°. 

Subject headings: Stars: Novae, Cataclysmic Variables, Stars:White Dwarfs, Stars: 
Individual: Constellation Name: QU Carinae 



Introduction 



Cataclysmic variables (CVs) are semi-detached binary stars in which a late-type m ain 
sequence star loses mass onto a white dwarf (WD) by Roche lobe overflow (jWarnerlll995l ). A 
sub-class, the NL systems (discussed below) are of interest because they are photometrically 
stable. It is believed that an analytic expression describes the radial temperature profile on 
the accretion disk, leading to a corresponding spectral energy distribution (SED). Tests using 
stellar SEDs (IWadelll988r) produ ced unsatisfactory results, but a theoretical model tailored 
to accretion disks (lHubenylll990l ) is available and might be expected to lead to satisfactory 
results. In this paper we show that synthetic SEDs, built on the theoretically-prescribed 
model, do not fit the observed data, but non-standard empirical models fit the observations 
closely. 

In non-magnetic systems the mass transfer stream produces an accretion disk with mass 
transport inward and angular momentum transport outward by viscous processes. The ac- 
cretion disk may extend inward to the WD; the outer boundary extends to a tidal cutoff limit 
imposed by the secondary star in the steady state case. If the mass transfer rate is below a cer- 
tain limit, the accretion disk is unstable and undergoes brigh tness cycles (outbursts), and if 
abov e the limit, the a c cretion disk is stable again st outbursts (IShafter. Wheeler, fc Cannizzo 
19861 : ICannizzd Il993l : IWarnen Il995l : lOsakil Il996l ). The latter case objects are called nova- 
like (NL) systems; their accretion disks are nearly fully ionized to their outer (tidal cutoff) 
boundary and this condition suppresses dwarf nova outbursts. NL systems are of special in- 
terest because they are expected to have an accretio n disk radial temperature profile given by 
an analytic expression (IFrank. King &: Raindll992l . eq.5.41)(hereafterFKR) thereby defining 
the so-called standard model. 



B ased on reports of QU Carin ae as an emission-line object (jStephenson. Sanduleak. fc Schild 

19681 ). iGilliland fc Phillips! (119821 ) (hereafter CP) undertook a spectroscopic study to deter- 



- 3 - 



mine orbital characteristics of the system. They determined an orbital period of 10.9 hours. 
No secondary star spectrum was detected, indicating that the spectrum is dominated by 
light from an accretion disk or the primary star. Balmer line emission is relatively weak. 
This evidence in turn suggested a high rate of mass transfer and a correspondingly bright 
system. It also suggests that QU Car is a NL-type system. We discuss known parameters 
of QU Car in more detail in §4. 

High speed photor netry confirni e d erratic varia t ions o f 0.1-0.2 mag. on time scales of 
minutes, as reported by ISchildl (Il969l ). iKnigge et al.l (Il994j ) used lUE spectra to study time 
variability of the UV resonance lines. They found that there are short-time-scale variations 
with a possible correlation with orbital period, but the data are insufficient to confirm the 
orbital nature of the variations. 



Hartley. Drew. &: Longl (120021 ) obtained Hubble Space Telescope STIS spectra of QU 
Car at three epochs. Blueshifted absorption in N V and C IV provides evidence of an 
outfiow wit h a ma ximum velocity of ~ 2000 km s~^. Further analysis of these spectra is in 
Drew et al.l (120031 ). The latter authors suggest that the QU Car secondary is a carbon star, 
and that QU Car is the most luminous NL system known, possibly at a distance of ?5i2000 
pc. 

Our objective in this study is to use synthetic spectra, based on a physical model of QU 
Carinae, to constrain properties of this interesting system. 



2. Interstellar reddening and the Hipparcos distance 



VerbuntI (119871 ) de termined E(B-V ) = 0.1±0.03 for QU Car from the 2200A interstellar 
medium(ISM) feature. iDrew et al.l (120031 ) used the STIS spectra (the same applied in this in- 
vestigation) to determine a neutral hydrogen column density of (6± 1) x 10^° cm~^. For a stan - 



dard gas-to-dust ratio of N{B.1)/Eb-v = 5.8x10^^ mag ^ cm ^ (IBohlin. Savage, fc Drake 
978[) . the corresponding reddening is E{B — V) = 0.1±0.015, in agreement with Verbunt 



(119871 ). In §8, we model the ISM atomic and molecular hydrogen absorption lines and find a 
hydrogen column density of A^(Hl) = 2.5 x 10^°cm~^, slightly lower than the value obtained 



by Drew et al., but within the range of acceptable values for a reddening of E(B-V)=0.1 (the 
Bohlin et al. empirical law has a very large scatter). We adopt E{B — V) = 0.1. 



The Hipparcos parallax (jPerryman et al.l Il997l ) of 1.64±1.50(mas), corresponds to a 
distance of 610 pc, with la errors giving a range between 318 pc and 7143 pc. Thus the 
Hipparcos parallax does not impose a strong constraint. 
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3. The FUSE, STIS, and lUE spectra 

The FU SE and STIS spectra were obtained at a time separation of about three years. 
A conversion factor was anticipated to match the two spectra, and we found it necessary to 
muhiply the STIS data by the factor 1.263 to match the flux level of the FUSE spectrum in 
the overlap region. Note also that the STIS spectrum used in our analysis is a combination 
of three exposures obtained over an interval of nearly one year. The data were dereddened 
to E(B-V)=0.1 (§2) after matching the flux levels. 

The FU SE spectrum of QU Car was obtained through the 30"x30" LWRS Large Square 
Aperture, in TIME TAG mode and was obtained on 23 April 2003, dataset 01560102000, 



PI N ame: L. Hartley. The data were processed with CalFUSE version 3.0.7 (jPixon et al. 



20071 ) which automatically handles event bursts. Event bursts are short periods during an 
exposure when high count rates are registered on one or more detectors. The bursts exhibit 
a complex pattern on the detector; their cause is as yet unknown (it has been confirmed that 
they are not detector effects). The main change from previous versions of CalFUSE is that 
now the data are maintained as a photon list (the intermediate data file - IDF) throughout 
the pipeline. Bad photons are flagged but not discarded, so the user can examine the filter 
and combine data without re-running the pipeline. 

The spectrum actually consists of 3 exposures, each covering a portion of a single orbit. 
These exposures, totaling 10,748 sec of good exposure time, were combined (using the IDF 
files). Consequently, the FUSE spectrum has a duration of about 0.27 of the period of QU 
Car. The spectral regions covered by the FUSE spectral channels overlap, and these overlap 
regions are then used to renormalize the spectra in the SiCl, LiF2, and SiC2 channels to the 
flux in the LiFl channel. We combined the individual channels to create a time-averaged 
spectrum with a linear 0.05A dispersion, weighting the flux in each output datum by the 
exposure time and sensitivity of the input exposure and channel of origin. In this manner 
we produced a final spectrum that covers the full Ft/^-E" wavelength range 905 — 1187A. 



The STIS spectra used in this study are those obtained in 2000 by lHartley. Drew, fc Long 



(120021). The three STIS spectra were processed with CALSTIS version 2.22. The H ST /ST IS 
spectra of QU Car were obtained in TIME TAG operation mode using the FUV MAMA 
detector. The observations were made th rough the 0.2x0.2 arcsec ^ aper ture with the E140M 



optical element. Additional details are in [Hartley. Drew, fc Longl (120021 ). The STIS spectra 



have a total duration of 7500 sec or about 0.2 of the period of QU Car. The STIS spec- 
trum, centered on the wavelength 1425A, consists of 42 echelle spectra that we extracted 
and combined to form a single spectrum extending from about 1160A to 1710A and with a 
resolution of O.lA. 
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lUE spectra extend the wavelength coverage to 3000A from the STIS hmit of 1710A, 
thereby providing a significantly enhanced constraint on a synthetic spectrum model. The 
lUE archive contains 6 SWP spectra and two LWP spectra of QU Car for the dates 26- 
27 June 1991, among others. Because of the intrinsic variability of the system, we have 
restricted our lUE data to that time interval. 

We chose a single pair of lUE spectra, SWP41927 and LWP20700, which were taken in 
immediate sequence on 27 June 1991. The close temporal sequence minimizes the possibility 
of source changes affecting the observations. We dereddened the spectra individually and 
merged them with the lUEMERGE facility, with the 'concatenate' option. We used an 
expanded ordinate scale to study the overlap of the STIS and merged lUE spectra and 
found that visual inspection provided a sensitive test of accurate superposition. In this way 
we found that the optimum divisor to superpose the merged spectrum on the scaled and 
dereddened STIS spectrum was 0.85. 

Figure 1 shows the FUSE and the STIS spectra of QU Car before dereddening. Note 
the very large number of interstellar lines in the FUSE spectrum. In Table 1 we identify the 
prominent absorption lines of the FUSE spectrum. The wavelength tabulations (column 
2) are the theoretical wavelengths. For clarity we did not include the multitude of ISM 
lines in t h e tab le. A complete list o f ISM lines (hydrogen + metals) can be found e.g. in 



Sembachl (Il999l ) and ISembach et al.l (l200ll ). Table 2 identifies the absorption lines in the 
STIS spectrum. The ISM lines of the STIS spectrum have been analyzed in Drew et 
al. (2003), and the broad absorption lines (ClllA1176, NvA1240, OvA1371, SilvA1398, 
ClvA1549, HellA1640) have been studied in detail for each STIS exposure separately in 
Hartley et al. (2002). The broad lines are seen in emission with a velocity broadening of 
several thousand km/s; an absorption component is superposed with a velocity broadening 
of a few hundred km/s. 



4. Parameters adopted initially 



The orbital period of 10.9 hours was determined by CP, using radial velocities from the 
AA4630 — 4660 emission feature and the A4686 Hell emission line. There is no detectable 



orbital light modulation (ISchildlll969l . CP), suggesting an inclination i<60°. The presence 
of doubled emission lines (CP) suggests i>30°. 



Based on the mai n sequence mass -radius relation of iLacyl (119771 ) and the usual period- 
mean density relation (IRobinsonl 119761 ) . CP obtained an approximate M2(= Ms) = 1.3M0. 
Note that this determination does not depend on an explicit knowledge of i. CP spectropho- 
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tometry sets a limit oi ^ 1% on the depth of any absorption feature from the secondary 
star in the QU Car optical spectrum. Comparison with the ~ 10% depth of absorption 
features in G-K standards (appropriate to the secondary star) establishes that the disk and 
primary of QU Car must be radiating at > IOLq. The strong He II A4686 emission, weaker 
He II A4542 emission, and the absence of He I lines at A4387, A4471, and A4713 suggests a 
temperature in the line-emitting region in excess of 20,000K. This is consistent with a fully 
ionized accretion disk expected for a NL system. GP determine a minimum distance to QU 
Car of 500 pc, based on = 11.2 for the system and sta ndard assumptions for interstellar 
absorption. From the My — Porb relation of I Warner! (jl995l . Figure 4.16), with My = 2.4 and 
interstellar absorption of 0.3 mag. (§2), we obtain a minimum distance that agrees with the 
GP value. 



A coarse estimate of the mass transfer rate is possible usin g the s tudy by ISmakI (119891 ). 



From Smak's Figure 2 we obtain M ~ 5 x lO~^M0yr~^. Later. ISmakI (119941 ) slightly revised 
the calculated mass transfer rates. This value is roughly c onsistent wi t h the mass transfer 
rate-orbital period relation of iPattersonI (Il984j . Figure 21). iDrew et al.l (120031 ) find that the 
QU Car C/He abundance may be as high as 0.06, an order of magnitude higher than the 
solar ratio, and the C/0 abundance is estimated to be greater than 1. Based on a power law 
fit to r nerged lUE spect ra, and a possible cutoff of the spectrum at 228A (the Hell series 
limit), brew et all J2003h suggest a mass transfer rate of M = 8.0x10 ^ MQyr ^ if QU Car 
is at a distance of 500 pc and M = 1.0x10^^ Moyr^^ if the distance is 2000 pc. 

GP cite the uncertainty in the Ms value and the lack of knowledge of i, and so do 
not derive a value of M^d- However, some constraint on M^d is possible as shown by the 
following discussion. 



From the standard expressions for radial velocity amplitude ( ISmartI 119471 . p. 360, eq. 
60), we obtain 

1 97rc:in7' 

(1) 



1 27rsinz ^ 



where we assume zero orbital eccentricity; the mass ratio q = Mg/My^d, i is the orbital 
inclination, K^^ is the WD radial velocity amplitude, P is the orbital period, and D is the 
component separation. From Kepler's third law. 



1 + 



47r2 



q GP^ M, 



(2) 



From these two equations we have 



P sinz 
' 27r K^d 



(3) 
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GP determine a value K^^ = 115 km sec~^ from He II emission lines (insecure because 
emission lines may not trace the WD motion accurately). GP derive an orbital period of 
0.4542'^. With these parameters determined by observation, an adopted value of i and an 
assigned value of determines D from equation (3), q from equation (1), and M^d from 
Mwd = Ms/q. Thus, a given i has an associated contour in the mass-mass plane that 
determines M„d values corresponding to assigned Mg values. 

Figure 2 illustrates the principle. The WD axis is bounded above by the Chandrasekhar 
mass. The diagonal hne represents the case of equal WD and secondary star masses; systems 
represented by points above this line are excluded by dynamical instability. (Mass transfer 
from a more-massive secondary causes the period to shorten, the separation of components 
and the secondary component Roche lobe to shrink, and so accelerates the mass transfer.) 
The value of i cannot be greater than about 60° to avoid eclipses. No models with i as 
small as 30° are permitted (for q = 1.0 the masses of both stellar components then are 
larger than the Chandrasekhar mass and for g < 1.0 the WD mass becomes even larger). 
A Ms > I.3OM0 is n ot permitte d because the secondary must lie on o r above the ZA MS 



eq. 



mass-radius relation (iLacyl 119771 ). The empirical mass-period relation of IWarnerl (1l995l . 
2.100) gives Ms = 1.29M0. A value Ms = 1.20Mq, with M^d = I.2OM0 is permitted with 
i = 38°. The uncertainty in leads us to adopt i = 40°, M^d = Ms = I.2OM0 as an 
appropriate case for simulation. 



5. Calculation of system models 



Our calculation of system models uses the BINSYN suite (ILinnell fc Hubenyl Il996 
Linnell et al.ll2007al ). The word 'model', as used in this section, has two context-dependent 
meanings. In one context it refers to a listing of the run of physical parameters in an an- 
nulus, or to properties of the collection of annuli comprising an entire accretion disk. In 
a different context it refers to the array of output synthetic spectra from BINSYN. In the 
latter context, orbital phase-dependent synthetic spectra separately of the system, the WD, 
the secondary star, the accretion disk, and the accretion disk rim, all produced in a single 
run, constitute a BINSYN model of a CV system. The synthetic spectra include effects of 
eclipse and irradiation. 



A BINSYN ac cretion disk rnodel uses an array of annulus models calculated wit 



program TLUSTY flHubenvlll988l . ll990l : lHubenv fc Lanj|l995l : iHubenv fc Hubenv 



1 the 



1998) and 



corre sponding synthetic spectra calculated with program SYNSPEC (IHubeny. Lanz. fc Jeffery 
19941 ). Table 3 is an illustration of an accretion disk model for M^d = I.2M0 and M = 
6.0x10"'' M0yr~^. The 26 TLUSTY annulus models cover the interval from the WD equa- 
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tor to the tidal cutoff boundary. The cutoff boundary is calculated from iWarnerl (119951 . p. 
57, eq. 2.61). Note that the annulus radii in column 1 are measured in units of a zero 
temperature WD. Successive columns provide some of the data available in the TLUSTY 
output for a given annulus model; each row describes a separate annulus model. The annuli 
all are H-He models, (i.e., these are the only explicit atoms; the next 28 elements are treated 
implicitly. See the TLUSTY manual for details. 0) Convection was ignored. The models 
converged for all of the Table 3 annuli. In some instances, for other M values, the outermost 
annuli failed to converge and we used so-called gray models (see the TLUSTY manual). 

The mo column lists the mass per unit area (gm cm~^) between the central plane and 
the upper boundary. Calculation of an annulus model specifies an initial uppermost layer 
in gm cm~^. Our calculation used 10~^ gm cm~^ for the initial layer. An annulus model 
typically may have 70 depth points. is the central plane temperature, log g is the value 
of that parameter at optical depth near 0.7, zh is the height, in cm, of the optical depth 0.7 
level above the central plane, and tross is the Rosseland optical depth at the central plane. 

SYNSPEC is used to calculate a synthetic spectrum for each TLUSTY annulus, with 
specified spectral resolution. The line list includes data for the first 30 periodic table ele- 
ments. Next, BINSYN is used with its specified number of annuli (different from the set of 
TLUSTY annuli) to calculate the various output synthetic spectra by integration over the 
array of BINSYN annuli for the accretion disk and by corresponding integrations over the 
stellar photospheric grids for the latter objects. Synthetic spectra for the individual BINSYN 
annuli follow by interpolation among the array of SYNSPEC spectra of the TLUSTY an- 
nuli. It is at this level that information is incorporated on orbital inclination, eclipse effects, 
radial velocity effects, etc., as needed for individual photospheric segments of all objects in 
the system. 

Calculation of a synthetic spectrum corresponding to a given M model is a computationally- 
intensive process and compromises are necessary in deciding how many TLUSTY annulus 
models to calculate to achieve acceptable accuracy in the final system synthetic spectrum. 
We believe our representational accuracy is high, but we cannot be sure that the fit in, e.g.. 
Figure 8 (following) could not be improved with a finer grid of annulus synthetic spectra for 
the TLUSTY annulus models and a larger number of BINSYN annuli. We note that this 
entire process must be repeated for each value of M and in turn for each value of M^d- 

In general, the WD synthetic spectrum contribution to the system synthetic spectrum 
is a strong function of the adopted WD T^q. Since the WD radius varies with the adopted 
WD Tefj, it makes a difference whether we use the zero temperature radius or the corrected 



^http:/ /nova. astro. umd.edu 
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radius in calculating annulus radii, in BINSYN, measured in units of the WD radius. We 
must use the corrected radius with BINSYN, designated by the symbol r^d, since the output 
spectra must describe the WD contribution for its Tcfj actually adopted. 

The WD Tes likely will be unknown in advance of the model calculation. We determine 
what annulus models to calculate with TLUSTY by using a fixed set of multiples of the 
zero temperature WD radius; we stress that, if we were to use the fixed multiples procedure 
with a variable WD radius, it would be nearly prohibitive to recalculate the entire array 
of TLUSTY annuli with each new approximation to an adopted WD Teg. Such iterative 
recalculation is unnecessary; BINSYN temperature-wise interpolation among the fixed (with 
'"wd = Twd,o) TLUSTY annulus models leaves the accretion disk synthetic spectrum essentially 
unaffected as the WD Teg- dependent radius changes. We adopt solar composition for the 
annulus models. In the present application the WD contribution is too small to matter, but 
it is important to retain the distinction, as this paper does. 



Our initial model, at all mass transfer rates, adopted a standard model (FKR. lLinnell et al 



2007al ). A standard model assumes the accretion disk is in hydrostatic equilibrium, that the 
accretion disk material is very nearly in Keplerian rotation, and that a single viscosity pa- 
rameter applies to t he entire accretion disk. A ll of our TLUSTY annulus models adopted a 



viscosity parameter (IShakura &: Sunyaevlll973l ) a = 0.1. 



6. QU Carinae models with M„d = I.2OM0, q = 1.00 

The long orbital period and the weak Balmer line emission, suggesting a NL-type sys- 
tem (GP), in turn indicate a hot WD. The radius o f a 1.2M(^ WD is slightly sens i tive t o 
the photospheric temperature. Table 12C-H-S.dat of iPanei. Althaus. fc Benvenuto J2000h . 



for homogeneous Hamada-Salpeter models, lists a radius r„d,o = 5.6 x 10 ^Rq for a zero 
te mperature, carbon, 1.2M(n WD- Ba sed on an adopted 55,000K WD, we used Figure 4a 



of iPanei. Althaus. &: Benvenuto! (I2OOOI ) to determine a corrected radius of 6.7xlO~^i?0. We 
will find that the WD makes an insignificant contribution to the system luminosity; the only 
impact the WD T^s has is indirect-temperatures in the accretion disk at specified multiples 
of the WD radius vary with the Tcff-dependent WD radius. 

A smaller assumed WD T^s produces a smaller i?„d and a correspondingly larger ac- 
cretion disk Teff in the innermost rings (FKR) since they then are in a deeper potential 
well. However, this effect is smaller than the temperature effect, on all annuli, of different 
assumed mass transfer rates. For uniformity in comparing different assumed mass transfer 
rates we arbitrarily adopt a fixed WD Teff=55,000K. The tidal cutoff radius of the accretion 
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disk is 150i?wd for the adopted WD Teg, and also equals O.998i?0. (Contrast with the same 
tidal cutoff radius, measured in units of the zero temperature WD radius, in the last line of 
Table 3). The BINSYN model uses 45 division circles on the accretion disk, with the first 
division at the WD equator, the next two at 1.1806i?wD5 and LSGll-R^d, and the last at the 
tidal cutoff boundary. 

We tested models with M = 0.1, 0.3, 0.6, 1, 3, 6, and 10x10"^ MqJt-^. In the first case a 
standard model produces an accretion disk temper ature profile whose outer annuli are cooler 



than 6000K. This model is unstable (jSmakI Il982l ) and must be rejected. In all remaining 



cases, the standard model produces an accretion dsik profile with a spectral gradient that is 
too large. 

Figure 3 illustrates the problem for M = 6.0x10"^ MqJy~'^] other values of M produce 
comparable plots after applying appropriate normalizing factors. In all cases in this paper, 
normalizing factors have been determined by eye estimates. In Figure 3 the peak flux of the 
WD, near lOOOA, is less than 0.01 of the system flux; the WD contribution to the system 
spectrum is negligible. Note that the synthetic spectrum consists of a continuum spectrum 
plus absorption lines of H and He II. (He II but not He I because He II is a single electron 
ion.) The spectral gradient problem persists if i is arbitrarily changed to 50 deg. For reasons 
discussed in §10, the M = 1.0x10"^ MqJt~^ rate is a maximum, so our range of M values 
cover the full range of possible choices for this M^d- 

We conclude that, for this WD and secondary star mass, no standard model SED 
provides an acceptable fit to the observational data. We believe the spacing between our M 
values is small enough to extrapolate our conclusion to any M value not actually calculated. 

Truncated standard models are of interest because the inner accretion disk can be elided 
by radiation from a hot WD or by interaction with a magnetic field associated with the WD. 
We ignore whether either process is likely for QU Car and restrict consideration to whether 
a truncated accretion disk can provide a SED that fits the observations with acceptable 
accuracy. 

A truncated accretion disk model, with the residual disk on the standard model temper- 
ature profile, should have a lower spectral gradient and so is a possible alternative. We have 
calculated truncated models for the last six mass transfer rates listed above (The model with 
the lowest rate of mass transfer failed because of its temperature behavior at large radii). 
Truncation cuts off the annuli that contribute most to the far UV. In the six largest M 
cases, the best truncation compromise either still had a too large spectral gradient or failed 
to represent the FUSE spectral region with acceptable fidelity. 

Figure 4 illustrates the fit for M = 6.0x10"'' MQyr~^, and a truncation radius of 
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-Rtrunc = 22.5i?wd- The overall fit is an improvement over Figure 3, but the fit to the FUSE 
spectrum is unsatisfactory and the spectral gradient still is too large. A larger truncation 
radius would improve the spectral gradient problem but make the fit to the FU SE spectrum 
worse. A smaller truncation radius has a still more discrepant spectral gradient. Other values 
of M are unsatisfactory for similar reasons. 

Having demonstrated that no standard model SED or truncated standard model SED 
fits the observational data, we consider non-standard model SEDs. The too-large spectral 
gradient can be reduced with a cooler (average) accretion disk. However, the flux maximum 
in the FUSE spectrum near lOOOA requires the presence of a hot region on the accretion 
disk. This suggests reducing the temperature of the disk model at intermediate radii. 

It is important to note that a modification of the temperature profile removes the direct 
connection to a specific mass transfer rate. Not only is the empirical profile disconnected 
from the standard model, the total flux from the empirical model differs from the flux of the 
parent standard model. We conducted very extensive preliminary experiments and found 
that a steep drop in the accretion disk temperature, followed by a much lower temperature 
gradient than in the standard model gave a fairly good fit to the observed spectra. However, 
we determined that further improvement was possible. 

It is desirable to preserve model connection to some mass transfer rate. On the standard 
model, 1/2 of the potential energy hberated by mass falling from the LI point (located on 
the secondary star photosphere) to the surface of the WD is converted to emitted radiation 
by the accretion disk (FKR). The remaining half either is liberated in a boundary layer, 
goes into spinning up the WD, or is deposited elsewhere, as in a wind. Based on energy 
considerations, a non-standard model accretion disk still can be associated with a mass 
transfer rate if roughly 1/2 of the total liberated potential energy appears as flux from the 
accretion disk. 

With that constraint in mind, we performed further tests. In the standard model the 
maximum accretion disk temperature occurs at -R(Tmax) = 1.3611i?wd- The region closer to 
the WD is the boundary layer region (discussed further in §10). Standard model temper- 
atures, for the M values we are considering, are so high at the temperature maxima that 
only the Rayleigh- Jeans tail contributes to the observed spectra from those regions. (If, in 
addition, as described above, 1/2 of the total hberated potential energy were deposited in 
the boundary layer region the Teflf there would become much higher.) The bolometric flux 
contribution, per unit area, of the innermost annuli far exceeds that of more remote parts 
of the accretion disk. Consequently we have preserved the standard model temperature 
maxima of the boundary layer region and have modified the temperatures of the remain- 
ing annuli. This choice preserves a connection with a particular M through the sensitivity 
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of the total accretion disk flux to the contribution of the highest temperature annuh. For 
convenience of nomenclature, we designate the region between the WD and the accretion 
disk R = 1.361 li?wd as the boundary layer even though this term has a slightly different 
connotation in the literature. 

It has proved possible to achieve excellent flts by adopting an isothermal proflle over 
almost all of the accretion disk while preserving a narrow region adjacent to the boundary 
layer, discussed in more detail below, and assigned an intermediate temperature value. The 
narrow region is important in achieving a fit to the short wavelength part of the FUSE 
spectrum. The highest temperature boundary layer region, important in the total accretion 
disk energy budget, contributes to the flux shortward of lOOOA in the observed spectrum by 
an amount that varies slowly with its temperature (the contribution is in the Rayleigh- Jeans 
tail) . Changing the boundary layer Teff from its standard model value produces appreciable 
changes in the system flux (thereby disconnecting from the adopted M designation) without 
producing a sensitive flt to the observed FUV flux. The narrow intermediate temperature 
region provides a reasonably sensitive means to produce a good flt to the observed FUV flux. 
In all of our models, the contribution of the intermediate temperature region near lOOOA also 
is in its Rayleigh- Jeans tail. 

The range of isothermal values that we have tested and that produce acceptable fits for 
this Mq, is hmited. The range extends from 15,000K at M = 1.0 xlO"'^ Mqyt-^ to 17,000K 
at M = 1.0x10"^ M0yr~^ (higher isothermal temperatures produce a too-large spectral 
gradient) . 

FigTuc 5 presents a sample fit to the observed spectra, again for M — 6.0x10"^ M0yr~^. 
The isothermal region has a temperature of 17,000K. In this comparison we have calculated 
a line spectrum with a resolution of 0.02A. The resolution of the FUSE spectrum is O.OSA. 
Our choice of synthetic spectrum resolution was to permit an accurate correction for ISM 
absorption, described in §8. Figure 6 shows details of the fit to the FUSE spectrum and 
the FUV end of the STIS spectrum. Extensive absorption lines of the ISM complicate the 
comparison. Further consideration of this topic is in §8. 

Figure 7 shows the fit to the STIS spectrum. The failure to fit Ly a, as we show 
subsequently, is rectified by including ISM absorption. The 1200A to 1500A interval is 
affected by the ISM. With the exception of not fitting the He II emission line, the fit in the 
1560A to 1710A interval, shown in Figure 8, generally is very good. Finally, Figure 9 shows 
the accurate fit to the red end of the SWP spectrum and the blue end of the LWP spectrum. 

The bolometric luminosity for the M — 6.0x10"^ Moyr~^ standard model accre- 
tion disk is 1.01x10^^ ergs sec~^. The bolometric luminosity of the empirical model is 
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3.95x10^^ ergs sec ^. Figure 5-Figure 9 are representative of corresponding plots for the 
four largest mass transfer rates for this WD and secondary star mass. 

The system parameters for this M^d are in Table 4. The default BINSYN model (be- 
fore imposing a non-standard accretion disk temperature profile) assigns standard model 
properties (FKR) at the various annulus radii, including the annulus semi-thickness. Those 
properties vary slightly with M but the tabulated value of within the thickness resolution 
listed, is appropriate for the range of M included in this study. Because our empirical models 
preserve a connection to corresponding standard model M values, the H value remains the 
same in the empirical models. 



7. QU Carinae models with M„d = O.GM©, q = 0.83 

We wish to test the possibility of acceptable standard model SED fits to the observed 
spectra over the full range of possible M^d and M values. Figure 2 indicates that the smallest 
possible Mwd for this system is about O.GMq. Again from Figure 2, M„d = O.GMp) and 



Ms = O.SMq fall on the i = 55° contour. From Figure 4a of iPanei. Althaus. fc Benvenuto 



fl2000h . for a M„d = O.GM© with T^s = 55, OOOK, we find i?wd = 1.59x10-27?©. The 
corresponding radius of a zero temperature WD is 1.21xlO^2-R0. For models with a range 
of M values, and the WD Teff=55,000K, the tidal cutoff radius of the accretion disk is 53R^d, 
which also is O.84Oi?0. 

We tested models with M = 0.3,0.6,1,3,6, and 10x10"^ Mgyr'^ The outer part of 
a standard model accretion disk for the first mass transfer rate is cooler than 6000K and is 



unstable ISmakI (119821 ). That model must be rejected. In all remaining cases, the standard 
model produces an accretion disk profile with a spectral gradient that is too large. The 
effect is similar to the M^d = I.2OM0 cases. Figure 10, illustrating the standard model for 
M = 3.0x10"'' MQyr~^, is representative of the standard model synthetic spectrum fits. As 
with the M = I.2OM0 case, we find that no standard model SED fits the observed spectra. 

Figure 11 presents a truncated model for M = 3.0x10^^ Moyr"^. The truncation 
radius is -Rtrunc = 9.91-Rwd- The overall fit is a big improvement over Figure 10 but the fit 
to the FUSE spectrum is unsatisfactory. The same type defect is present for the other M 
values. We conclude that no truncated model provides a SED which fits the observed spectra 
satisfactorily. As with the Mwd = I.2OM0 case, it proves possible to produce good quality 
synthetic SED fits for the four largest mass transfer rates we tested. Figure 12 illustrates 
the result for M = 3.0x10^^ M0yr~^. Fits for other mass transfer rates are of similar 
quahty. Figure 13 presents the FUSE fit, and Figure 14 the STIS fit. Figure 15 shows 
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the expanded scale fit, 1560A to 1710A. Figure 16 shows the fit to the lUE spectra and the 
accurate connection to the red end of the STIS spectrum. Parameters of the model for this 
M^d are in Table 5. The comments describing Table 4 also apply to Table 5. 

In summary, over the full range of permissible M^d masses in this system, no acceptable 
M has an associated standard model SED that fits the observed spectra with acceptable 
accuracy. The same conclusion follows for truncated standard models. As with the M^d = 
1.2Mq case, all permissible M cases have an associated empirical accretion disk model that 
provides an accurate fit to the set of observed spectra. Further discussion is in §10 



8. Inclusion of a model of the ISM 



We have developed software to simulate the ISM and have used this routine to determine 
ISM properties in the hne of sight to QU Car. A custom spectral fitting package is used 
to estimate the temperature and density of the interstellar absorption lines of atomic and 
molecular hydrogen. The ISM model assumes that the temperature, bulk velocity, and 
turbulent velocity of the medium are the same for all atomic and molecular species, whereas 
the densities of atomic and molecular hydrogen, and the ratios of deuterium to hydrogen 
and metals (ii icluding helium) to hydrogen are adjustab le parameters. The model use s 
atomic data of iMortonI (120001 . l2003l ) and molecular data of lAbgrall. Roeuff. fc Drira J2000h. 
The rn olecular hydrogen transmission values have been checked against those of iMcCandliss 
(120031 ). Most of the lines in the FUSE spectrum are due to molecular hydrogen. 



The ISM pa r amete rs are in Table 6. The HI column density is in reasonable agreement 
with iDrew et al.l (120031 ). Although Drew et al. present evidence for a two-component ISM 
toward QU Car, we have restricted our simulation to a single component. The output of 
the simulation program is a list of transmission values at the same spectral resolution as the 
system synthetic spectrum. The product of the ISM transmission function and our O.O2A- 
resolution synthetic spectrum produces a corrected synthetic spectrum. We have used the 
M = 6.0x10"^ M0yr~^, M^d = I.2OM0 case to illustrate inclusion of the ISM simulation. 



The ISM-corrected synthetic spectrum line density is so great in the FUSE region 
that it is not meaningful to present a plot comparable to Figure 6 (and a color plot is not 
helpful). Figure 17 shows a partial plot of the FU SE comparison; note the cutoff at 912A as 
compared with Figure 6. The ISM lines in our new model fit the observed lines well in many 
cases, including line depths. There are depth discrepancies in some cases, and our model 
misses some observed lines. Our ISM model is a reasonable first order representation of a 
very complex ISM. There now is a good fit at Ly 7 and Ly 6. Large unmodeled residuals 
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remain from N IV and S VI (see Figure 1). iHartley. Drew. &: Long (120021 ) discuss features 
like these in terms of a wind and/or a disk chromosphere. No adjustment of parameters in 
our present set of models can produce a fit to these features, and they represent the largest 
remaining residuals. Note the differences in the ordinates in Figure 17 and Figure 1. The 
latter apply to unreddened fluxes and the former to fluxes corrected for reddening. 

Figure 18 shows the fit near Ly /?. Note the unmodeled broad O VI features at 1033A and 
1037A and the improved Ly 13 fit as compared with Figure 6. Figure 19 shows the fit to 
the ST IS spectrum. The fit to Ly a now is accurate (compare with Figure 7). The ST IS 
spectrum shows a number of unmodeled absorption features including the N V line pair near 
1243A and the series of absorption lines starting near 1370A. 

Although we have illustrated correction for the ISM for a single empirical model, the 
procedure is applicable to all of those models. Our ISM model fit to the C I line s in the 
FU SE spectrum does not indicate a C overabundance in the ISM. iDrew et al.l (120031 ) present 
evidence of a C overabundance in the matter undergoing transfer from the secondary, based 
on C emission lines. 



9. Determination of system parameters 



Several conclusions follow from the analysis to this point. 

(1) No standard model that we studied fits the observed spectra in either of the two 
Mwd cases. For reasons discussed previously, we believe the same conclusion would be true 
for any intermediate case, and we believe the two cases are limiting cases. 

(2) No truncated standard model that we studied produces an acceptable fit to observed 
spectra in either of the M^d cases. Again, for reasons already discussed, we believe that 
situation would be true for any intermediate case, and that the two cases are limiting cases. 

(3) Possible M values are bounded for small values by an unstable accretion disk. This 
limit is about M = l.OxlO"'^ MqJi'^ 

(4) A possible constraint at high M is t heoretical: expansion of the WD to a red giant 
and c onsequent engulfing of the entire binary (jPaczyhski fc Zytkowlll978l : ISion. Acierno. fc Tomczyk 
19791 ). Ho wever, the studies cited co nsidered spherical accretion and do not apply to disk 



accretion (jPaczyhski fc Zvtkow 



1978). The effects of mass transfer onto WDs via accretion 



disks have been studied by lPiro fc BildstenI (120041 ). based on a spreading layer model. These 
authors consider M values up to 1.5xlO~^MQyr~^ and find no significant expansion of the 
WD. Note that their study considers the mass actually transferred to the WD, and in the case 
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of QU Car there is a wind which may carry a way an appreciable mass ( iHartley. Drew, fc Long 



20021). At M = 1.0x10"^ Meyr"^ and larger IPiro fc BildstenI (120041 ) point out that hydrogen 



is expected to be burning steadily on the WD surface, leading to luminosities in the range 
10^^ — lO^^erg s~^. A luminosity that large exceeds the disk luminosity and so would be in- 
compatible with the observed double A4686 emission line. We take M = 1.0x10"^ M^yT'^ 
to be an upper limit to acceptable mass transfer rates. 

(5) Non-standard models can be found that produce acceptable fits to the combined 
observed spectra. We achieved successful fits with a hot accretion disk region bracketed 
between the WD and a radius R = 1.3611i?wd (which we designate as a boundary layer), 
an intermediate temperature region adjacent to the hot region, and an isothermal region ex- 
tending to the tidal cutoff boundary. The hot region Teg was set equal to the standard model 
temperature maximum to preserve a designation connection with the standard model M. 
Because the peak emissions of both the hot region and the intermediate temperature region 
are at much shorter wavelengths than the short wavelength limit of the FUSE spectrum, 
we have no good way to impose tight temperature constraints on those regions. 

Table 7 compares properties of models that fit the observed spectra well. Columns 3,4,5 
list the Teff values for the boundary layer, the intermediate region, and the isothermal region. 
The last four lines list the possible range of isothermal region Tes values. The first four lines, 
for Mwd = I.2M0, present, first, a comparison of line 1 with line 6 for the same Tes values 
of the intermediate and isothermal regions. Line 6, for M„d = O.6M0, corresponds to a 
slightly smaller accretion disk (compare in Table 4 and Table 5) and the value of i for 
Table 4 is 40° versus 55° for Table 5. Line 3 and line 4 preserve the same value of isothermal 
region Tes and illustrate the increase of bolometric luminosity with, primarily, the change 
in the boundary layer Teg. Note that the calculated distance to QU Car is the same for 
these two lines, indicating that the normalizing factor to fit the synthetic spectrum to the 
observed spectra depends strongly on the isothermal region Teg and is fairly insensitive to 
the boundary layer T^s. Column 6 lists the bolometric luminosity of the standard model 
for comparison with the bolometric luminosity of the empirical model in column 7. Note 
that the entries in the two columns are roughly comparable, in agreement with our earlier 
discussion on preserving a connection between the empirical model and the standard model 
M. The last column lists the distance to QU Car for each empirical model. The range of 
distances is not large for M^d = I.2M0, and all are in rough accord with the Hipparcos 
distance (§2). The smallest distance for M^d = O.6M0 (332 pc) is in marginal disagreement 
with the Hipparcos distance. 

There is no clear distinction, among the acceptable models, of the quality of fit. It is 
important to note that the FUSE, STIS and lUE spectra were obtained at different times. 
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that they were fitte d together by visual estimate, an d that the QU Car luminosity varies on 
a short time scale (ISchildl Il969l : iKnigge et al.lll994l ). Because of this continuing variation, 
a new study with new data likely would find larger differences from the present data than 
the residuals in our present model fits. Although accurate fits in some spectral regions are 
possible (e.g., Figure 8), the fit accuracy should not be overinterpreted. 

The fit in the 920A to 1200A region is fairly sensitive to the Teg of the intermediate 
temperature region. The width of that region arbitrarily was set to two (BINSYN) annuli, 
and their corresponding widths depend on the choice of 45 annulus divisions of the accretion 
disk in the BINSYN model. By eye estimate, the fits with the present models are optimal. 
We suspect that a choice of one of the models over another, based on, say, a test, could 
arise from numerical effects (interpolation approximations, etc.) rather than from a closer 
representation of physical reality. Since the largest remaining residuals arise from unmodeled 
wind or chromospheric features, these models present the closest representation of the QU 
Car data currently available. 

Table 4 and Table 5 list the orbital parameters used for the two values of M^d- Since 
there is no light curve to simulate, the Teg of the secondary star is of importance only to 
the extent that it could affect the long wavelength end of the system synthetic spectrum. A 
nominal 6500K secondary in the M„d = I.2M0 system provides no evidence of contamination 
in the system synthetic spectrum. Initial retention of that Tcfj in the M^d = O.6M0 system, 
followed by substitution of a 3500K secondary produces no detectable effect on the system 
synthetic spectrum. 

We are confident that M^d masses could be selected between our limiting values of 
O.6M0 and I.2M0, combined with appropriate values of Ms and i from Figure 2, assigned 
a range of M values and produce SED fits to the observed spectra comparable in quality 
to Figure 5 or Figure 12. Table 7 represents the extremes of possible system models, based 
on existi ng spectrosco pic data. Because of its consistency with the empirical mass-period 
relation (IWarnerlll995l . eq. 2.100) we favor the M^d = I.2M0 solution. 



10. Discussion 



We have appropriated the term "boundary layer" to designate the region between the 
WD and the accretion disk R = 1.3611i?wd- This region more generally is associated with a 
very hot region occupying the sa me space but predicted to emit hard X-rays if optically thin 
or soft X-rays if optically thick (j Warned 1 19951. sect. 2.5.4). This prediction is modified if a 



wind carries off part of the accretion mass (IHoare fc Drewlll991 



I993I ). In the latter case. 
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predicted boundary layer Teg's are no higher than the values adopted in our models. Based 
both on its geometric location and the much higher local temperature than in the rest of 
the accretion disk, our appropriated term actually is in keeping with the usual sense of the 
term. 

Adoption of the intermediate temperature layer obviously is a crude approximation. 
A temperature discontinuity at its inner and outer boundary, as in our present models, 
is not physically credible. However, this model is a much better fit to observation than 
with our earlier experimental (§6) models, which adopted a sharp, nonetheless more gradual 
temperature drop to the isothermal disk region than in this paper. 

The outermost accretion disk region contibutes little to the system synthetic spectrum. 
The Tefj values, for a limited region which would vary from one M to another, could be closer 
to standard model values than the isothermal assumption without producing a detectable 
change in the system synthetic spectrum. On the other hand it is known that temperatures 
in the outer region of an accret ion disk need an upward correcti on from the standard mo del 
to account for impact heating (lLasotall200ll ) and tidal heating (IBuat-Menard et al.ll200ll ). 



The accretion disk temperature profile s found in this study differ from profiles found 



for oth er NL systems. In both SDSSJ0809 flLinnell et al.ll2007af ) and IX Vel flLinnell et al. 



2007bl ) the final model included an inner isoth ermal region with a n outer region that followed 
a standard model. In the case of MV Lyr (ILinnell et al.l |2005| ) in a high state, the final 
model had a standard model accretion disk extending from an inner truncation radius to 
an intermediate radius with an isothermal radius beyond. It is of interest that all of the 
NL systems included in these studies have accretion disks whose temperature profiles clearly 
differ from the standard model. 



11. Summary 

The analysis in this study leads to the following conclusions: 

(1) Standard model SEDs, within the tested ranges of M and M„d, are too hot to fit 
the combined set of FU SE, STIS, and lUE spectra of QU Car with acceptable accuracy. 

(2) Truncated standard model SEDs fail to fit the combined set of FU SE, STIS, and 
lU E spectra of QU Car for similar reasons. 

(3) Alternative non-standard models have SEDs that accurately fit the set of observed 
spectra. These models individually have a boundary layer T^s equal to the maximum ac- 
cretion disk Tefj of an associated standard model M. A narrow adjustable T^s annulus is 
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needed adjacent to the boundary layer to "fine tune" the SED fit to the FUSE spectrum. 
The remainder of the accretion disk is isothermal to within the accuracy of the SED fit. The 
model-to-model residuals arc so similar that it is not possible to choose a best model. The 
QU Car WD mass is greater than 0.6Mq and less than the Chandrasekhar limit; a solution 
with M„d = I.2M0 is consistent with an empirical mass-period relation. The M from the 
secondary star is between l.OxlO~''M0 yr~^ and l.OxlO~^M0 yr~^. The set of models agree 
on a limited range of possible isothermal region Teff values between 14,000K and 18,000K. 

The orbital inclination is between i = 40° and i = 60°. The Hipparcos distance, 610 pc, 
is representative of the empirical model results. 

(4) A model of the ISM combines with the empirical accretion disk models to provide 
an improved fit to the observed FUSE and STIS spectra. The largest remaining residuals 
are from unmodeled, broad, high excitation features probably associated with a wind or a 
disk chromosphere. 

(5) The accretion disk dominates the system flux in the FUSE, STIS, and lUE spec- 
tral regions. An assumed 55,000K WD provides a negligible flux contribution, as does the 
accretion disk rim and the secondary star. 

The authors are grateful to the referee for a careful and prompt consideration of the 
initial version of this paper. The final version has benefitted substantially from the criti- 
cisms of the referee. PG is thankful to Mario Livio for his kind hospitality at the Space 
Telescope Science Institute. Support for this work was provided by NASA through grant 
number HST-AR-10657.01-A to Villanova University (P. Godon) from the Space Telescope 
Science Institute, which is operated by the Association of Universities for Research in As- 
tronomy, Incorporated, under NASA contact NAS5-26555. PS is supported by HST grant 
GO-09724.06A. 

This research was partly based on observations made with the NASA/ESA Hubble 
Space Telescope, obtained at the Space Telescope Science Institute, which is operated by the 
Association of Universities for Research in Astronomy, Inc. under NASA contract NAS5- 
26555, and the NASA-CNES-CSA Far Ultraviolet Explorer, which is operated for NASA by 
the Johns Hopkins University under NASA contract NAS5-32985. 

REFERENCES 

Abgrall, H., Roueff, E., & Drira, I. 2000, A&AS, 141, 297 
Bohlin, R. C., Savage, B. D., & Drake, J. F. 1978, ApJ, 224, 132 



-20- 



Buat-Menard, V., Hameury, J.-M., & Lasota, J.-R, 2001, A&A, 366, 612 
Cannizzo, J. K. 1993, ApJ, 419, 318 
Dixon, W. v., et al, 2007, PASP, 119, 527 

Drew, J. E., Hartley, L. E., Long, K. S., & van der Walt, J. 2003, MNRAS, 338, 401 

Frank, J., King, A. & Raine, D. 1992, Accretion Power in Astrophysics (Cambridge: Univ. 
Press) (FKR) 

Gilliland, R. L. & Phillips, M. M. 1982, ApJ, 261, 617 

Hartley L. E., Drew, J. E., & Long, K. S. 2002, MNRAS, 336, 808 

Hoare, M. G. & Drew, J. E., 1991, MNRAS, 249, 452 

Hoare, M. G. & Drew, J. E., 1993, MNRAS, 260, 647 

Hubeny, I. 1988, Comp. Phys. Comm., 52, 103 

Hubeny, I. 1990, ApJ, 351, 632 

Hubeny, I., Lanz, T., & Jeffery, C. S. 1994, in Newsletter on Analysis of Astronomical Spectra 
No. 20, ed. C. S. Jeffery (CCP7;St. Andrews: St. Andrews Univ.), 30 

Hubeny, I., & Lanz, T. 1995, ApJ, 439, 875 

Hubeny, I., & Hubeny V. 1998, ApJ, 505, 558 

Knigge, C, Drew, J. E., Hoare, M. G., & la Dous, C. 1994, MNRAS, 269, 891 

Lacy, C. H. 1977, ApJS, 34, 479 

Lasota, J.-P. 2001, New Astronomy Reviews, 45, 449 

Linnell, A. P., & Hubeny, I. 1996, ApJ, 471, 958 

Linnell, A. P., Szkody, P., Gansicke, B., Long, K. S., Sion, E. M., Hoard, D. W., & Hubeny, 
I. 2005, ApJ, 624, 923 

Linnell, A. P., Hoard, D. W., Szkody, P., Long, K. S. Hubeny, I., Gansicke, B., & Sion, E. M. 
2007a, ApJ, 654, 1036 

Linnell, A. P., Godon, P., Hubeny, I., Sion, E. M., & Szkody P. 2007b, ApJ, 662, 1204 



- 21 - 



McCandliss, S. R. 2003, PASP, 115, 651 
Morton, D. C. 2000, ApJS, 130, 403 
Morton, D. C. 2003, ApJS, 149, 205 
Osaki, Y. 1996, PASP, 108, 39 

Paczynski, B., & Zytkow, A. N. 1978, ApJ, 222, 604 

Panei, J. A., Althaus, L. G., & Benvenuto, O. G. 2000, A&A, 353, 970 

Patterson, J. 1984, ApJS, 54, 443 

Perryman,M. A. C. et.al. 1997, A&A, 323L, 49P 

Piro, A. M., & Bildsten, L. 2004, ApJ, 610, 977 

Robinson, E. L. 1976, Ann. Rev. Astr. Ap., 14, 119 

Schild, R. E. 1969, ApJ, 157, 709 

Sembach, K. R., in Stromlo Workshop on High Velocity Clouds, ASP Conference Series, 
Vol.666, 1999, B. K. Gibson and M. E. Putman (eds) ( |astro-ph/9811314^ available 
from the FUSE website) 

Sembach, K. R., Howk, J. C, Savage, B. D., Shull, J. M., & Oegerle, W. R., 2001, ApJ, 561, 
573 

Shafter, A. W., Wheeler, J. C, & Cannizzo, J. K. 1986, ApJ, 305, 261 

Shakura, N. I., & Sunyaev, R. A. 1973, A&A, 24, 337 

Sion, E. M., Acierno, M. J., & Tomczyk, S. 1979, ApJ, 230, 832 

Smak, J. 1982, Acta Astron., 32, 199 

Smak, J. 1989, Acta Astron., 39, 317 

Smak, J. 1994, Acta Astron., 44, 45 

Smart, W. M. 1947, Spherical Astronomy (Cambridge: University Press) 
Stephenson, C. B., Sanduleak, N., & Schild, R. E. 1968, Ap. Letters, 1, 247 
Verbunt, F. 1987, A&AS, 71, 339 



-22- 

Wade, R.A. 1988, ApJ, 335, 394 

Warner, B. 1995, Cataclysmic Variable Stars (Cambridge: University Press) 



This preprint was prepared with the A AS IM^jX macros v5.2. 



-23- 
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Siv 1073.1+1073.5 


s 


2.0 


unresolved doublet 


1081.85 


Fell 1081.87 


c,ism 




affected by H2 



-24- 



Table 1 — Continued 



Wavelength 

(A) 


Line 
Identification 


Origin 


FWHM 

(A) 


Comments 


1083.97 


Nil 1083.99 


c,ism 


0.15 




1096.85 


Fell 1096.88 


c,ism 


0.1 


affected by H2 


1112.00 


Fell 1112.03 


c,ism 


0.1 




1118.05 


Siv, VI 1117.9+Pv 1117.98 


s 


1.7 




1121.97 


Fell 1121.97 


c,ism 


0.1 




1122.40 


Si IV 1122.48 


s 


1.5 




1 1 oo on 


r V iizo.Ui+oiIV iizo.od 


s 


i.o 




1133.70 


Fell 1133.67 


c,ism 


0.1 




1134.15 


Nl 1134.17 


a,c,ism 


0.1 




1134.45 


Nl 1134.42 


a,c,ism 


0.1 




1135.00 


Nl 1134.98 


a,c,ism 


0.15 




1142.40 


Fell 1142.37 


c,ism 


0.1 




1143.25 


Fell 1143.23 


c,ism 


0.1 




1144.95 


Fell 1144.94 


c,ism 


0.1 




1152.85 


Pll 1152.82 


c,ism 


0.1 




1168.75 


Hel 1168.9 


a?s? 


0.3 


emission 


1172.50 


Si 1172.55? 


c,ism 


0.3 




1173.00 


Oiii 1173.37 


s? 


0.3 




1175.95 


Cm 1174.9-1176.4 


s 


1.7 


shifted by +0.35A 



Note. — a=possibly including air glow contamination; c=circumbinary material; 
ism=ISM; s=source 
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Table 2. STIS Lines 



Wavelength 


Line 


Origin 


FWHM 


Comments 


(A) 


Identification 




(A) 






N T 1 1 68 54 


*^,lk5lll 


1 

U. X 




1 1 71 OR 


Nt 1171 OR 


IGTTl 
L-^lblll 


n 1 ^ 

U. X(J 




1 1 75 7 


Cttt 

Vv' 111 


q 

ij 




imT'pcnlAT'pri mnl'I'Tnlp'l' 

Lllll CoWl Vvyvl 111U.1U1 L/ICL 




n T 1 1 88 83 


L/.lolll 


W. X 






Si TT 11 90 90 

kJl 11 XXiyW.Z.W 


V/.lOlll 


07 

W. W 1 




1 1 Qf) 9Q 


n T 1 1 90 95 

1 XXiyW.ii-W 


L/. lolll 




n n T* 1 at'^^H 

U.111 CoVJl V CLl 


X X ijyj .Tnj 


Si TT 11 90 49 

kJl 11 XXtyW.T:.^ 


L/.lOlll 


13 

W. X^J 




X X ijij .yjij 


C T 1 1 93 03 

1 XXt^tJ.WO 


^^.lOlll 


07 

W . W 1 






Si TT 1 193 29 

kJl 11 XXt^O.^t^ 


*^.ioiii 


16 

W. X w 




1 1 Q7 91 


Mn TT 1 1 97 1 8 

iVXll 11 J. J-Cf 1 . XO 


*^,lk5lll 


05 




1 1 qq 4"? 


Mn TT 1 1 99 39 

IVXll 11 X X (J (_/ 


'^.lOlll 


05 

W. W(J 




1 1 QQ 5fi 


N T 1 1 99 55 


'^.lOlll 


16 

W. X w 




1 9nn 9'^ 


Nt 1900 99 

ly 1 x.^ww..^.^ 


*^,lk5lll 


14 

W. Xrt 




1 9nn 79 

X ^ WW . 1 ^ 


Nt 1900 71 

IN 1 X^WW. 1 X 


L.. lolll 


1 3 

W. ±kJ 




1201.15 


Mn TT 1901 19 

IVXll 11 x^wx.x^ 


c isixi 


05 

w • wu 




1 906 5S 

X WW . tJt_> 


Si TTT 1 906 50 

kJllll XZjWW.WW 


L.. lOlll 


19 

w . X^ 




1915 7 

XZj J.0 . 1 


Ht 191 5 66 

XXI XZ(XtJ.WW 


o,L^,lolll 


19 




1 9'^7 nq 


npTT 1 937 06 


*^,lkMll 


06 






N V 1 938 89 




9 9 




1 9'^q qfi 


Mp-tt 1939 93 


*^,1&111 


05 




1/4U.42 


Mgll 1240.39 


c,ism 


0.05 




1243.00 


Nv 1242.80 


s 


1.7 




1250.62 


Sll 1250.58 


c,ism 


0.11 




1253.83 


Sll 1253.81 


c,ism 


0.12 




1259.53 


Sll 1259.52 


c,ism 


0.12 




1260.45 


Sill 1260.42? 


? 


0.26 


asymmetric 


1260.77 


Ci 1260.74 


c,ism 


0.06 




1266.4 


Cl 1266.42 


s 


1.6 


shallow 


1272.7 


On 1272.23 


s 


1.5 


shallow 


1277.29 


Znll 1277.31? 


c,ism 


0.08 
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Table 2 — Continued 



ITT 1 J 1 

Wavelength 


Line 


(Jrigm 


FWHM 


Comments 


(A) 


T 1 1 ' C 1 ' 

Identiiication 




(A) 






C, T 1 977 95-1-1 977 98? 


*^,lk5lll 


08 




1 977 CC7? 


Ct 1977 51-1-1977 55? 


*^,lolll 


19 

yj. j-Zj 


vciy oiictiivjw, U.111 Coiji V CLi 




n T 1 ocf) 14 


*^,lolll 


05 


q n f» 1 1 f^wT 

ollctllij w 




P TT 1 301 87 

± 11 X Ov^ X . O ( 


'^.lolll 


05 




X OW^ . X 


O T 1 309 1 7 


*^,lolll 


90 




1 304 SQ 


Qi TT 1 '?n4 37 


L/. lolil 


14 

\J . Xt: 




1 sn4 8^ 

X OWt:. (JtJ 


O T 1 304 86 


L.. lolll 


90 


vci y oiidiiwvv 


1 '^Dfi 09 


O T 1 306 03 


L^.lOlil 


20 


V Cl V olidjllW VV 




Ni TT 1 31 7 99 


*^jlolll 


09 


oiidinj w 


1 "^98 87 




*^jlk5lll 


06 


q n 1 1 r^AUT" 

kMldllVJ w 


1 '^9Q 14 

X kJ^U • Xrt 


Ct 139Q 19 


'^.lolll 


06 


AT'PTAT' *3 n Pi 1 1 mxT" 
V CI y oiicuiu w 


1334 56 


C II 1334 53 


c ism 


0.21 




1 335 73 


n TT 1 335 71 

11 X OOO . I X 


*^,lolll 


14 

yj. Xt: 




1 347 98 


CI T 1 347 24 

V^l 1 XOt: 1 


V^.lOlil 


05 


c n fi 11 n^AT" 

ollfXllU VV 


1 355 65 


O T 1 355 60 


k^.lOlil 


05 


\rov\r t;rii^llow 
V y oiicxiiwvv 


1 358 81 

X O^JO . (J X 


Cii TT 1 358 77 

Li 11 X O^JO . 1 1 


L/. lolil 


05 


\rfiT\r Qri allow 
V y oiidiiiv^vv 


1 37n 1 7 


Ni TT 1 "^70 14 

IM 11 XOiU.X^ 


L^.lolll 


07 


c n d 1 1 rwnT 
olicLllU W 


1 371 6 

X (J 1 X . w 


O V 1 371 30 


q 

ij 


1 6 

X . W 




1 3Q3 Q5 


Si TV 1 3Q3 76 

kJl IV A-kJUkJ. 1 \J 


q 


1 1 

X . X 




1 403 Q5 


Si TV 1 409 77 

kJl 1 V ±rt\J^. 1 1 


c 

ij 


1 9 




1 /I on T 


O 

1 


s: 


1.7 


seen only in one exposure 


1414.43 


Gall 1414.40 


c,ism 


0.05 


very shallow 


1454.89 


Nill 1454.84 


c,ism 


0.04 


very shallow 


1526.73 


Sin 1526.71 


c,ism 


0.17 




1532.58 


Pll 1532.53 


c,ism 


0.05 


very shallow 


1548.18 


Civ 1548.19 


s 


1.0 




1550.85 


Civ 1550.77 


s 


1.0 




1560.36 


Ci 1560.31 


c,ism 


0.06 




1560.74 


Ci 1560.71 


c,ism 


0.07 


very shallow 


1640.52 


Hen 1640.47 


s 


1.0 
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Tablc 2 — Continued 



Wavelength 


Line 


Origin 


FWHM 


Comments 


(A) 


Identification 




(A) 




1656.98 


Ci 1656.93 


c,ism 


0.07 




1657.43 


Cl 1657.38 


c,ism 


0.04 


very shallow 


1657.96 


Cl 1657.91 


c,ism 


0.04 


very shallow 


1670.81 


All 1670.79 


c,ism 


0.13 





Note. — a=possibly including air glow contamination; 
c=circumbinary material; ism=ISM; s=source 
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Table 3. Properties of accretion disk with mass transfer rate M = 6.0x10 ^ MqJt ^ and 

WD mass of 1.20 Mr,,. 



r/rwd,o 


en 








c 


log g 


Zu 


^Ross 


1. 


,36 


334084 


2, 


.213E5 


4. 


,534E6 


7, 


.88 


6.92E7 


8. 


,73E4 


2. 


,00 


299540 


2, 


.701E5 


4. 


,221E6 


7, 


.62 


1.21E8 


1. 


,02E5 


3, 


,00 


242217 


2, 


.640E5 


3. 


,386E6 


7, 


.31 


1.99E8 


9. 


,93E4 


4. 


,00 


203586 


2, 


.458E5 


2. 


,802E6 


7, 


.08 


2.76E8 


9. 


,34E4 


5. 


,00 


176588 


2, 


.284E5 


2. 


,394E6 


6, 


.89 


3.56E8 


8. 


,80E4 


6. 


,00 


156665 


2, 


.133E5 


2. 


,097E6 


6, 


.75 


4.38E8 


8. 


,36E4 


7. 


,00 


141312 


2, 


.004E5 


1. 


,871E6 


6, 


.62 


5.20E8 


7. 


,99E4 


8. 


,00 


129082 


1 


.894E5 


1. 


,693E6 


6, 


.51 


6.04E8 


7. 


,68E4 


9. 


,00 


119082 


1 


.798E5 


1. 


,550E6 


6, 


.42 


6.90E8 


7. 


,43E4 


10. 


,00 


110733 


1 


.714E5 


1. 


,432E6 


6, 


.33 


7.77E8 


7. 


,22E4 


12. 


,00 


97540 


1, 


.572E5 


1. 


,247E6 


6, 


.18 


9.54E8 


6. 


,88E4 


14. 


,00 


87537 


1, 


.460E5 


1. 


,110E6 


6, 


.06 


1.14E9 


6. 


,63E4 


16. 


,00 


79657 


1, 


,366E5 


1. 


,004E6 


5, 


.95 


1.32E9 


6. 


,45E4 


18. 


,00 


73268 


1, 


.288E5 


9. 


,197E5 


5, 


.85 


1.51E9 


6. 


,31E4 


20. 


,00 


67967 


1, 


.221E5 


8. 


,505E5 


5, 


.77 


1.71E9 


6. 


,21E4 


24. 


,00 


59649 


1 


.110E5 


7. 


,439E5 


5, 


.62 


2.11E9 


6. 


,08E4 


30. 


,00 


53387 


1 


.024E5 


6. 


,334E5 


5, 


.50 


2.73E9 


6. 


,02E4 


35. 


,00 


45436 


9, 


.088E4 


5. 


,682E5 


5, 


.32 


3.27E9 


6. 


,03E4 


45. 


,00 


37854 


7, 


.930E4 


4. 


,782E5 


5, 


.12 


4.37E9 


6. 


,20E4 


60. 


,00 


30685 


6, 


.767E4 


3. 


,957E5 


4, 


.89 


6.11E9 


6. 


,61E4 


80, 


,00 


24852 


5, 


.762E4 


3. 


,304E5 


4, 


.66 


8.53E9 


7. 


,35E4 


100. 


,00 


21092 


5, 


.081E4 


2. 


,891E5 


4, 


.49 


I.IOEIO 


8. 


,21E4 


120. 


,00 


18440 


4, 


.581E4 


2. 


,601E5 


4, 


.34 


1.37E10 


9. 


,15E4 


140. 


,00 


16458 


4, 


.195E4 


2. 


,383E5 


4, 


.22 


1.64E10 


1. 


,02E5 


160. 


,00 


14911 


3 


.885E4 


2. 


,211E5 


4, 


.11 


1.91E10 


1. 


,12E5 


180. 


,00 


13667 


3 


.631E4 


2. 


,069E5 


4, 


.01 


2.19E10 


1. 


,22E5 
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Note. — Each lin e in th e table represents a separate annulus. A 
Shakura fc SunyaevI (119731 ) viscosity parameter a = 0.1 was used in 



calculating all annuli. r^d,o is the radius of a zero temperature, car- 
bon, Hamada-Salpeter WD. See the text for additional comments. 
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Table 4. QU Car Model System Parameters, M^d = 1.2M^ 



parameter 


value 


parameter 


value 




I.2OM0 (adopted) 


rs(pole) 


1.185i?0 


M, 


1.2Mq 


(point) 


1.664i?0 


P 


0.4542 days 


rs(side) 


1.245i?0 


D 


3.3275i?0 


r<iback) 


1.348i?0 




500.0 


log (pole) 


4.39 




3.75 


log (point) 


-5.00 


i 


40° (adopted) 


log Qsiside) 


4.30 




55, OOOK(nommal) 


log (back) 


4.15 




O.OO666i?0 


ra 


O.998i?0 


log g^d 


8.9 


n 


O.OO666i?0 






H 


O.O5i?0 



Note. — wd refers to the WD; s refers to the secondary 
star. D is the component separation of centers, is a Roche 
potential. specifies the outer radius of the accretion disk, 
set at the tidal cut-off radius, and is the accretion disk inner 
radius. H is the semi-height of the accretion disk rim (standard 
model) . 
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Table 5. QU Car Model System Parameters, M^d = 0.60M( 



parameter 


value 


parameter 


value 




O.eOM© (adopted) 


rs(pole) 


O.8748i?0 


M, 


O.SOM© 


Ts (point) 


1.235i?0 


P 


0.4542 days 


rs(side) 


O.9169i?0 


D 


2.56556i?0 


r^back) 


O.9977i?0 




162.2 


log fifs(pole) 


4.27 




3.47364 


log 5's (point) 


0.43 


i 


55° (adopted) 


log gsiside) 


4.19 


TeS,wd 


55, OOOK(nommal) 


log g's(back) 


4.03 




O.O159i?0 


ra 


O.84i?0 


log g^d 


7.8 


n 


O.O159i?0 






H 


O.O5i?0 


Note. — 


wd refers to the WD; 


s refers to the secondary 



star. D is the component separation of centers, is a Roche 
potential. Va specifies the outer radius of the accretion disk, 
set at the tidal cut-off radius, and r;, is the accretion disk 
inner radius. H is the semi-height of the accretion disk rim 
(standard model). 
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Table 6. Parameters of ISM model 



Parameter Value 



N(H) 2.5xl020cm-2 

D/H 1x10-^ 

N(H2) 7xl0i^cm-2 

vel. 0.0 km 

temp. 350 K 

metallicity 0.01 (1.0 — solar) 
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Table 7. Properties, empirical mass transfer models 



M 


M 




TeS 




Tcs 


-^bol 


Lhoi 


D 


(Me) 


(Mo yr-i 


) 


(b.l.,K) 


(inter., K) 


(isoth.,K) 


(std.) 


(empir.) 


(pc.) 


1.20 


1.0x10- 


-7 


185,000 


90,000 


15,000 


1.68x10^6 


1.06x10^^ 


550 


1.20 


3.0x10" 


-7 


240,000 


100,000 


16,000 


5.05x10^^ 


2.18x10=^^ 


610 


1.20 


6.0x10- 


-7 


290,000 


110,000 


17,000 


1.01 xlO^^ 


3.95x10^6 


687 


1.20 


1.0x10- 


-6 


330,000 


110,000 


17,000 


1.68x10''^^ 


5.52x10^^ 


687 


0.60 


1.0x10 


-7 


100,355 


70,000 


14,000 


2.56x10^^ 


2.77x10^5 


332 


0.60 


3.0x10- 


-7 


132,000 


90,000 


15,000 


7.68x10^^ 


7.11x10-^5 


406 


0.60 


6.0x10 


-7 


157,000 


100,000 


17,000 


1.54x10^^ 


1.22x10-^'^ 


475 


0.60 


1.0x10 


-6 


178,000 


110,000 


18,000 


2.56 xlO^*^ 


1.86x10^^ 


518 



Note. — b.l.= boundary layer; inter. = intermediate region; isoth.= isothermal region; 
std.= standard model; empir. = empirical model The units of columns 6 and 7 are erg s~^. 
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1 650 



1700 



Fig. 1. — F^/^ii^ plus 5'T/5' combined spectra with line identifications. The ^T/^ spectrum 
has been multiplied by 1.263 to match the flux level of the FUSE spectrum in the overlap 
region. 
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WD mass(M0) 

Fig. 2. — Mass-mass plot showing regions where a system model is permitted. The region 
above the diagonal line is excluded by dynamical instabihty. The region below i — 60° is 
excluded by the absence of observed echpses. See the text for details. 
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"TTT ^ 




1000 1500 2000 2500 3000 

Wavelength (A) 

Fig. 3. — Standard model fit to FUSE + STIS + lUE spectra. The synthetic spectrum is 
a continuum spectrum + H and He II lines. The mass transfer rate is M = 6.0x10"^ MqJI'^ 
and the WD mass is 1.2QMq. The synthetic spectrum spectral gradient is too large. Note 
that the ordinate scale is logarithmic. 
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1000 1500 2000 2500 3000 

Wavelength (A) 

Fig. 4. — Truncated model fit to observed spectra. The mass transfer rate is M = 
6.0x10"'^ MQyr"-"^, the WD mass is 1.20Mq, and the accretion disk has been truncated 
at Rtrun = 22.5i?wd- The spectral gradient is a better fit than Figure 3 but still unacceptable 
and the synthetic spectrum fails to fit the FUV accurately. 
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1000 1500 2000 2500 3000 

Wavelength (A) 

Fig. 5. — Comparison to observed spectra of model with 17,000K accretion disk face and 
associated M = 6.0x10^^ M^yi-K The WD mass is I.2OM0. The synthetic (line) spectrum 
(heavy line) has a resolution of O.O2A. The synthetic spectrum has been divided by 4.5x10^^. 
The corresponding distance to QU Car is 687pc. See the text for details. 
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900 1000 1100 1200 1300 

Wavelength (A) 

Fig. 6. — Detail of Figure 5. The synthetic spectrum (heavy line) fits the peaks of the 
FU SE spectrum (grey plot) but the latter spectrum shows strong effects of ISM absorption. 
The fit includes the short wavelength end of the ST IS spectrum (light line). See the text 
for details. 
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1200 1300 1400 1500 1600 1700 

Wavelength (A) 

Fig. 7. — Detail of Figure 5 showing fit to tfie STIS spectrum. ISM absorption affects tfie 
ST IS spectrum shortward of 1500A. The lUE contribution has been suppressed to avoid 
confusion with the STIS spectrum. See the text for a discussion. 
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5.0 r 




2.0 L-i 1 I I I I I I I I I I I I I . II 1 1 I I lilLi I 1 U I I I ] 

1560 1580 1600 1620 1640 1660 1680 1700 

Wavelength (A) 

Fig. 8. — Expanded scale detail of Figure 5, synthetic spectrum resolution 0.02A. The grey 
plot is a segment of the ST IS spectrum, resolution 0.1 A. Note that the ordinate scale is 
linear. Except for a few details and the He II emission line, the fit is quite accurate. Note 
the doubling of the He II line. The synthetic spectrum normalizing factor is the same as in 
Figure 5 (i.e., the normalizing factor has not been changed to produce a better fit in this 
spectral region). 
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1700 



1750 



1850 



1900 



1800 
Wavelength (A) 

Fig. 9. — Detail of Figure 5 showing fit of the synthetic spectrum (heavy hne) to the long 
wavelength end of the SWP lUE spectrum (hght slowly undulating hne) and the short 
wavelength end of the LWP lUE spectrum (light jagged line). The accurate superposition 
of the SWP spectrum on the red end of the ST IS spectrum (short, heavy noisy line) is 
visible at the extreme left of the plot. 



1000 1500 2000 ^ 2500 3000 

Wavelength (A) 

Fig. 10. — Fit of standard model continuum synthetic spectrum to combined observed 
spectra. The mass transfer rate is M = 3.0x10"^ MqJt:''^ and the WD mass is 0.6Mq. The 
synthetic spectrum spectral gradient is too large. Compare with Figure 3. 



-44- 




1000 1500 2000 2500 3000 

Wavelength (A) 

Fig. 11. — Fit of truncated model continuum synthetic spectrum, WD mass of O.QMq, 
to combined observed spectra. The mass transfer rate is M = 3.0x10-"^ MqYt-^. The 
truncation radius is -Rtrunc = 9.91i?wd- The synthetic spectrum fails to fit the observed 
spectra in the FUV. Compare with Figure 4. 
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1000 1500 2000 2500 3000 

Wavelength (A) 

Fig. 12. — Fit of an isothermal model synthetic spectrum (heavy line) to combined observed 
spectra. The synthetic spectrum has a resolution of 0.02A. The mass transfer rate is M = 
3.0x10"^ M0yr~^ and the WD mass is O.GMq. The synthetic spectrum has been divided by 
1.57x10^^. The corresponding distance to QU Car is 406 pc. See the text for details. 
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Fig. 13. — Detail of Figure 12 showing comparison with the FUSE spectrum (grey plot) 
and the short wavelength end of the STIS spectrum. Note the strong ISM absorption hues 
in the FUSE spectrum. 
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Fig. 14. — Detail of Figure 12 showing comparison with the ST IS spectrum. ISM absorption 
affects the ST IS spectrum shortward of 1500A. The lUE contribution has been suppressed 
to avoid confusion with the STIS spectrum. 
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Fig. 15. — Expanded scale detail of Figure 12. Note that the ordinate scale is hnear. 
The grey hne is a segment of the STIS spectrum, resolution O.lA. The synthetic spectrum 
resolution is 0.02A. Except for a few details and the He II emission line, the fit is quite 
accurate. Compare with Figure 8. The synthetic spectrum normalizing factor is the same 
as in Figure 12. 
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Fig. 16. — Detail of Figure 12 showing the fit of the synthetic spectrum (heavy line) to the 
long wavelength end of the SWP lU E spectrum (hght slowly undulating hne) and the short 
wavelength end of the LWP lUE spectrum (light jagged line). The accurate superposition 
of the SWP spectrum on the red end of the STIS spectrum (short, heavy noisy line) is 
visible at the extreme left of the plot. 
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Fig. 17. — ISM-corrected version of Figure 6. The spectral resolution of the synthetic 
spectrum is 0.02A. Strong absorption hues of N IV and S VI (grey underlying plot) remain 
unmodeled, possibly from a wind or accretion disk chromosphere. See Figure 1 for identifi- 
cations. There now is a good fit to Ly 7 and Ly S. Note the hnear ordinate scale. See the 
text for a discussion. 
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Fig. 18. — Comparison of ISM-corrected synthetic spectrum and FUSE spectrum near 
Ly 13. The ISM correction provides a good fit to a number of the narrow ISM hues. Note 
the hnear ordinate scale. Note the remaining unrepresented broad O VI hues at 1033A and 
IO37A. These features probably associate with a wind or chromosphere not included in the 
model. 
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Fig. 19. — Synthetic spectrum fit, including an ISM correction, to STIS spectrum. The 
heavy line is the synthetic spectrum, the grey line is the observed spectrum. The lUE 
contribution has been suppressed to avoid confusion with the STIS spectrum. There now 
is a good fit to Ly a, in contrast to Figure 7, but unmodeled features remain. See the text 
for a discussion. Note the linear ordinate scale. 



